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We study the link between relaxation to the equilibrium and anomalous superdiffusive motion in 
a classical N-body hamiltonian system with long-range interaction showing a second-order phase- 
transition in the canonical ensemble. Anomalous diffusion is observed only in a transient out-of- 
equilibrium regime and for a small range of energy, below the critical one. Superdiffusion is due 
to Levy walks of single particles and is checked independently through the second moment of the 
distribution, power spectra, trapping and walking time probabilities. Diffusion becomes normal at 
equilibrium, after a relaxation time which diverges with N. 
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In the recent years there has been an increasing inter- 
est for physical phenomena which violate the central limit 
theorem such as anomalous diffusion and Levy walks. 
These violations are not an exception in Nature and have 
been observed in many different fields and also in connec- 
tion with deterministic chaos in low dimensional systems 
[Q-^. The availability of more powerful computers has 
made possible to study deterministic chaos and subdif- 
fusive motion in systems with many degrees of freedom 
using nearest-neighbour coupled symplectic maps |^ . In 
a very recent work superdiffusive motion has been found 
in a N-body Hamiltonian system with long-range cou- 
plings 0. The mechanism underlying this anomalous 
diffusion is similar to the one proposed by Geisel in 
"egg-crate" two-dimensional potentials. 

In this Letter we present a novel study of superdiffusion 
and Levy walks in a Hamiltonian system of N fully cou- 
pled rotors (called Hamiltonian Mean Field, HMF) [|,||. 
The new interesting result is that, in HMF superdiffu- 
sion is connected to the presence of quasi-stationary non- 
equilibrium states, rather than to the mechanism pro- 
posed by Geisel ||l| and found also in Q . HMF has been 
used to investigate relaxation to thermodynaniical equi- 
librium for systems with long-range interactions. It has 
been studied both at a macroscopic level, by means of 
the canonical formalism and at a microscopic dynam- 



ical level. The canonical ensemble predicts a second- 
order phase transition from a clustered phase to a ho- 
mogeneous one 1^,^. On the other hand, microcanoni- 
cal simulations show a strong chaotic behavior in the re- 
gion below the critical energy; Lyapunov Exponents and 
Kolmogorov-Sinai entropy reach a maximum at the crit- 
ical point Q . These results have been confirmed also for 
long but finite-range interactions pO| . Of particular im- 
portance for this Letter are the results obtained in Ref. |^ 
concerning the discrepancies between microcanonical re- 
sults and canonical predictions. In fact, numerical simu- 
lations performed at constant energy reveal the existence 
of out-of-equilibrium Quasi-Stationary States (QSS) with 
an extremely slow relaxation to equilibrium. In Ref. |Tl| ] 
these QSS are shown to become stationary solutions in 
the continuum limit. 

The main results of this letter are: 

1) we find evidence of an anomalous superdiffusive be- 
havior below the critical energy. Anomalous diffusion 
changes to normal one after a crossover time Tc, as also 
found by other authors 0J|,^,0,^2|jl3j . Power spectra 
confirm the presence of the anomaly; 

2) the superdiffusive behavior is connected to the pres- 
ence of out-of-equilibrium QSS . We give substantial nu- 
merical evidence that the crossover time Tc coincides with 
Tr, the time needed for QSS to relax to canonical equi- 
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librium; 

3) we give an interpretation of our results in terms of 
Levy walks, that are originated by chaotic transport of 
each rotor, which moves with an energy not constant in 
time and alternately sticks to the cluster, which has a 
quasi-regular motion or undergoes free walks far from it 
with a constant velocity much greater than that of the 
cluster. Trapping time and walking time probability dis- 
tributions show a power law behavior. The corresponding 
exponents can be related to the superdiffusion exponent 
using the model of Ref. Q and are very similar to those 
found in the fluid flow experiment of Solomon et al. |^ . 

In the following we remind the formalism and then 
we review the numerical results. HMF describes a sys- 
tem of N classical particles (or rotors) characterized by 
the angles 9i and the conjugate momenta pi. Each ro- 
tor interact with all the others according to the following 
Hamiltonian: 

H{0,p)^K + V, (1) 

where 

i—l — 1 

are the kinetic and potential energy. One can de- 
fine a spin vector associated to each rotor nij = 
[cos{9i)T siniOi)] and a total magnetization M = 
-ir^^-^mi. The Hamiltonian then describes N classi- 
cal spins similarly to the XY model. This system has a 
ferromagnetic second-order phase transition from a clus- 
tered phase to a homogeneous one at a critical tem- 
perature Tc = 0.5 and a corresponding critical energy 
= Ec/N = 0.75 (see ref. |,§). The equations of 
motion for the N rotors are given by: 

n ^ Pi, Pi ^ -Msin{6i - (f)) , i = l,...,N , (3) 

where (M, (jj) are respectively modulus and phase of the 
total magnetization vector M. These equations are for- 
mally equivalent to those of a perturbed pendulum. To 
study relaxation to canonical equilibrium, we solve these 
equations on the computer using fourth order symplectic 
algorithms (the details can be found in ref. ^). We start 
the system in a given initial distribution and we com- 
pute 6i,pi at each time step, and from them the total 
magnetization M and temperature T (through the rela- 
tion T = 2 < K > /N). We consider systems with an 
increasing size N and different energies U = E/N . 

Diffusion and transport of a particle in a medium or 
in a fluid flow are characterized by the average square 
displacement cr^(t) in the long-time limit. In general, 
one has 

a^{t)^e (4) 



with a ~ 1 for normal diffusion. All the processes with 
a ^ I are termed anomalous diffusion, namely subdiffu- 
sion for < a < 1 and superdiffusion for I < a < 2. 

In order to study anomalous diffusion in HMF we fol- 
low the dynamics of N rotors starting the system in a 
"water bag" , i.e. a far-off-equilibrium initial condition 
obtained by putting all the rotors at = and giving 
them a uniform distribution of momenta with a finite 
width centered around zero. We compute the variance of 
the one-particle angle 9 according to the expression 

aUt) {e- <9>f> , (5) 

where < . > indicate the average over the N particles, 
and we fit the value of the exponent a in eq. (^. In fig.l 
we plot on a log-log scale ag vs. t for N — 500 at three 
different energies: U=0.2, U=0.6, U=5. The continuous 
lines are shifted fits and show a very clear power law over 
a few decades; the corresponding values for the slope a 
are indicated in the inset. The numerical results show 
clearly three different types of behavior: 

1) No diffusion for very low energy, i.e. U < 0.2. In 
this case all the particles belong to a single cluster and 
Q 0. 

2) A ballistic regime a = 2 for J7 bigger than the crit- 
ical energy {Uc = 0.75) (a short-time ballistic regime is 
obviously always present for all energies). 

3) Superdiffusion with a = 1.38 ± 0.05 for U=0.6, in 
the transient regime. After a crossover time ~ 7-10* a 
change to the slope a = 1 (normal diffusion) is observed. 
The superdiffusive regime is present in the energy range 
0.5 <U < 0.75. 

In fig. 2 we study the dependence on N of the anoma- 
lous diffusion and the coincidence of crossover time Tc 
with the relaxation time r^, i.e. the time the system 
needs to reach the canonical temperature (horizontal dot- 
ted line in panel b)). We report ag and temperature 
vs. time for = 500,2000 and U = 0.69. A slope 
a — 1.42 ± 0.05 is observed in a first time stage, in which 
the temperature is different from the canonical value. In 
fact the temperature maintains for a very long period 
a constant value which corresponds to a QSS belonging 
to the continuation of the homogeneous phase at a sub- 
critical energy (see in particular fig.l of Refs. and 
Ref. [|ll|). Indeed, the crossover time from anomalous 
to normal diffusion Tc coincides with the relaxation time 
Tr. This result has also been checked, changing the accu- 
racy of the numerical simulation. The transient regime, 
in which QSS and anomalous diffusion are present, in- 
creases linearly with N ||], consequently for N=2000 one 
gets superdiffusion over almost 3 decades. On the other 
hand, the slope a does not seem to strongly depend on N 
and moreover in the range 0.6 < U < 0.69 it varies from 
1.38 to 1.42. A similar scenario has been recently conjec- 
tured by Tsallis for systems with long-term memory 
and slow relaxation to equilibrium, but to our knowledge 
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this is the first time that it has been found in a numerical 
simulation. 

The importance of noise and finite-size fluctuations in 
the crossover from anomalous to normal diffusion has 
been studied in detail in Rcfs. [ p2| , p"3| j7|] . On the contrary, 
Kaneko and Konishi ||^ claim that relaxation to normal 
diffusion is due to phase-space uniform sampling, which 
occurs asimptotically. Our results show that this relax- 
ation occurs in coincidence with relaxation to equilibrium 
of QSS, which is a quite close mechanism to the one pro- 
posed in Ref . |6| . However, at variance with these latter 
authors, our model displays superdiffusion, rather than 
subdiffusion, in the transient. Some important physical 
facts might be crucial in the observed differences. Su- 
perdiffusion occurs near a second-order phase transition 
in our case, while in the model of Ref. ^ no phase tran- 
sition is present. The model of Ref. |Q has a flrst-order 
phase transition and the particles performing correlated 
flights belong only to a distinct dynamical class or phase 
for N oo. In our case this is not true: close to the 
critical energy fluctuations are maximal and do not dis- 
appear in the thermodynamical limit. Each particle reg- 
ularly performs free walks and trapped oscillations until 
it forgets the initial condition and tends to a brownian 
motion. 

Evidence in favour of this mechanism is provided by 
the link of superdiffusion with Levy walks. For low- 
dimensional chaotic Hamiltonian systems, superdiffusion 
has been interpreted as due to the trapping of the par- 
ticles by the cantori of the phase space; particles can 
eventually escape and walk freely before a new trapping 
occurs and this mechanism prevents from normal 
diffusion. An analogous situation occurs in the experi- 
ment of Solomon et al. for chaotic transport in a two- 
dimensional rotating flow In this case tracer parti- 
cles are trapped and untrapped by a chain of six vortices. 
This last mechanism is very similar to ours. In fig. 3 we 
report the time behavior of the angle 6 and of the cor- 
responding conjugate momentum p of a test particle in 
the transient anomalous diffusion regime (panels (a) and 
(b)) and in the equilibrium regime (panels (c) and (d)) 
for U — 0.6 and N = 500. Free walks and trapped mo- 
tion are observed in the transient regime; the walks have 
an almost constant velocity corresponding to the sepa- 
ratrix between bounded and free motion (~ of 
the perturbed pendulum of Eqs. (^. In the equilibrium 
regime the test particle remains trapped in the cluster; 
it oscillates around its center and drifts together with it 
on a much longer time scale (for a study of cluster mo- 
tion see Q). It is important to notice that the energy 
of the test particle is not conserved. The particle walks 
freely when accidentally it receives enough energy that 
allows it to escape from the mean field. In this sense the 
mechanism of anomalous diffusion in our case is similar 
to that of non-conservative systems. A quantitative dif- 
ference between the two behaviors can be obtained by 



performing the power spectrum of the motion in fig. 3(a) 
and 3 (c). We get a power law with slope —2 for the equi- 
librium regime, as it should be for brownian motion, and 
a slope smaller than —2 for the transient. To study the 
connection between Levy walks and anomalous diffusion 
we evaluate trapping and walking time distributions. A 
free walk is identified by A6' > 2n. In fig. 4 we consider 
for N — 500 and two energies U=0.6 and U=0.69, the 
probability distribution of "walking times" and "trap- 
ping times" . They show, as expected, a clear power law 
decay: 

Pwalk^t-l", Ptrap-t-" . (6) 

The values of and v obtained from the fitting are re- 
ported in figure. Their value is crucial because, the two 
exponents /i and v can be related to the anomalous dif- 
fusion coefficient a. The following relationships, derived 
in ref. Q, are the most appropriate for HMF: 

a^2 + v~ 11 2</i<3, j/<2 (7) 
a = 4-^ 2<^<3, i^>2. (8) 

These formulae are valid for a one-dimensional system 
under the assumptions of walks with a constant velocity, 
separated by sticking events with no motion. In the case 
shown we get for U=0.69 (0.6) = 2.14 ±0.1 (1.98 ± 
0.1), V = 1.58 ± 0.05 (1.34 ± 0.05) and thus a value for 
a = 2+v-ii = 1.44±0.05 (1.36±0.05) which is consistent 
with the values obtained from the fits shown in figs.l and 
2, within our numerical accuracy (the relative error on 
the exponent a obtained by fitting the slope of the vari- 
ance is ~ 5%). As found in ref. ||^ we get for the trapping 
probabilities the exponent v < 2, which is not the usual 
case encountered in low dimensional conservative maps. 
The reason of that is likely the non-conservation of energy 
for the test particle motion. In conclusion, we have found 
superdiffusion and Levy walks in a Hamiltonian system 
showing a second-order phase transition. This behavior 
occurs in a transient out-of-equilibrium regime in a range 
of energies sligthly smaller than the critical energy where 
the system is strongly chaotic and QQS exist. We have 
also found that the equilibration time to reach the canon- 
ical temperature, which diverges with N, corresponds to 
the crossover to normal diffusion confirming a recently 
proposed scenario ||l^. This feature has been observed 
for the first time in a deterministic chaotic system with 
many degrees of freedom and could be of relevance to un- 
derstand more realistic situations such as the anomalous 
diffusion observed in fluid flow experiments |^ . 
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FIG. 1. For N=500 we show three different behaviors; 1) 
no diffusion for U=0.2, 2) ballistic diffusion for U=5 and 3) 
superdiffusion for U=0.6. In this last case we considered the 
average over 5 events. The straight full lines are shifted fits 
and the relative slopes are also indicated. The relative errors 
obtained from the fits are ^ 5%. We show in the inset the 
numerical evaluation of the slope a vs time for the case U=0.6. 
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FIG. 2. Variance and temperature for two different sizes 
N=500 (dashed lines), 2000 (fuU lines) at U=0.69. Panel (a) 
shows that a ~ 1.4 does not depend on N (within the ac- 
curacy of the calculations) and occurs only in a transient 
regime. Once the canonical temperature, shown in panel (b), 
is reached, diffusion becomes normal. The relaxation time 
Tr is clearly larger for bigger A^. The vertical dashed lines 
indicate Tc ~ Tr for N=500 and N=2000. 
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FIG. 3. Angle 6 and momentum p of a typical particle for 
U = 0.6 and A'' = 500. Panels (a) and (b) refer to the tran- 
sient regime, (c) and (d) to the canonical equilibrium state. 
See text for more details. 




FIG. 4. For the cases U=0.6 and U=0.69, wc show the 
probability distribution functions for trapping (open dia- 
monds) and walking times (open circles) calculated in the 
transient regime (2000 < t < 8000). The correspondent fits 
and exponents are also indicated, see text for more details. 
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